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Abstract 

Energy spectra of precipitating electrons are fitted to the suo of three 
distributions, a power lav, a Maxwellian and a Gaussian. This fitting pro» 
ccdure detemines seven parameters idilch characterize the essential features 
of each spectrum. These characteristic parameters are used to carry- out 
various studies involving precipitating electrons. It is shown that the absence 
of the power 'lav population from a particular spectrum is related to the soft* 
ness of the precipitating primary flux, that the Maxwellian temperature and the 
Gaussian peak energy have a positive correlation the strength of which varies 
with local time, that the upward moving Gaussian population has a loss cone 
distribution, and that the one-dimensional velocity distribution parallel to the 
magnetic field occasionally displays a plateau or a hump on the tail. 



1. Introduction 




The precipitation of charged particles at high latitudes represents one 
of the rtost interesting phenomena occurring in Che ionosphere, because of the 
many interactions which these particles have undergone and which influence the 
shape of their spectra. A careful study of the characteristic features of their 
spectra is of special importance because those features represent the signature 
wf their region of origin and of their acceleration mechanism and therefore 
contain infotmation about physical processes occurring large distances away from 
the point of observation, frequently even deep in the magnetosphere. In addition, 
the precipitation spectrum determines the effect of the particles on the local 
ionosphere, such as ionization, emission of light, heating and the excitation 
of plasma instabilities. 

Since the advent of in-situ measurements oy instruments onboard space 
vehicles probably hundreds of millions of precipitating electron spectra have 
been obtained and duly stored on tape, ^^uch valuable information contained in 
those spectra rejoins to be extracted. One practical way to deal vith such 
a large data base is to associate with each spectrum a limited number of charac- 
teristic parameters and to base the study of the physics of precipitating elec- 
trons on those parameters. Obviously the choice of such charact eristic parameters 
is critical for they must contain the important physical Infomation carried by 
each spectrum. The parameterization of the electron fluxes lends itself not only 
to the interpretation of already existing data but can also be usefully included 
in the reduction of data currently being collected. The parameters of the various 
populations tell the investigator immediately the characteristic features of the 
observed spectra. 

It is the purpose of this paper to describe a method for obteining such a set 
of parameters for each spectrum and then to report the results of some studies. 



The basis for this method is the reeogaltioa that the observed fluxes 
generally represent the superposition of several different populations and 
that each of these populations can be described in tanas of certain character* 
istie parameters of an analytic function vbich is fitted to the observed fliaces. 

ft 

This has been done before by other investigators (e.g. Frank and Acker*son» 

1971; Burch et al., 1976; Lin and Hoffman, 1979]. To our knovledge, hcvever, 
no systeaatic effort has so far been made to parameterize large numbers of 
spectra or. a production basis by means of a computerized fitting procedure 
and to use the parameters for investigations of the physical processes under* 
gene by the flimses. 

In Section II the analytic expression selected by us is discussed. The 
characteristic parameters of esch flux spectrum are determined by the fit of 
said analytic function to the measured spectrum. The computerized fitting 
procedure is briefly described in Appendix: I. In Section III some cf the 
present studies are briefly described. Section IV contains a summary of 


the results. 




II. Analytic Representation of Sleetron Energy Spectra 

It has been pointed out in Section I that in order for a paraaetric repre* 
sentatien of precipitation fluxes to be useful it is necessary that the selected 
paraneters contain the important physical inforcation carried by each spectrum. 
Based partly on earlier work by other authors and after extensive trials involving 
various different models ve came to the conclusion that the electron spectra 
investigated by us so far can best be represented as a superposition of three 
populations given by the following expression: 
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Eere C ■ (2m ) is a known constant and a, c, n, T , A, E , and A are 
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the parameters determined by the particular shape of each spectrum. The first tens 
represents the well-known power-law spectrum which generally describes the combined 
backscattered and secondary electrons. Most of the electrons in the energy 
range from 200eV (the lowest energy measured by the Low Energy Electron Analyzer 
onboard the AE-D satellite) to IkeV frequently fall into this category. The 
second term is a Maxwellian with temperature T, while the third term 
represents a non -Maxwellian peak of Gaussian shape. Such peaks are characteristic 
of inverted-V spectra. The Gaussian distribution which has not been used in 
earlier works was istroduded because the Maxwellian distribution gave only poor 
agreement with spectra containing near-monoenerget ic peaks. If a Maxwellian fit 
were used for such spectra, the error was consisteritly large.' than for spectra 
without such peaks. Tne details of the fitting procedure are presented in 
Appendix I. The analytic expression given by equation (1) is not necessarily 
the optimum model for all energy ranges, altitudes or even instruments. Vt are 
continuing to test different models (for example, to allow for two maxima ) . 



The three populations which make up the resultant energy spectrum are believed 
to represent different physical processes undergone by the respective electrons. 

For example, the Gaussian tern may be the result of a magnet ie»f laid ^aligned 
electric field accelerating the electrons into a near-monoenergetic beam. It 1s 
therefore logical to investigate the properties of such fields in teztts of the 
parameters associated with the Gaussian rather than in terms of the total flux. 
Similarly, other investigations are best carried out by studying the properties 
of the individual populations, including of course any correlations between them. 

Figures la and lb show two examples both ef which have Gaussian peaks. 

Yet, as these exa^les indicate, each of the three populations of equation (1) 
is not necessarily represented in every spectrum. Tne reason for the presence 
(or absence) of any of the three populations may shec light on the origin of 
the precipitating electrons and also on the Interactions they have undergone. 
Figure la is a typical auroral spectrum, while the st-ectrus of Figure lb almost 
looks like £ magnetoshe?ith spectrum, both as regards its shape as veil as the 
value of the peak energy. Since the altitude of the spacecraft was at 610 km, 

Che spacecraft was obviously not in the magnetosheath. However, Fester and 
Burrows [1977] have recently proposed chat magnetosheath plasma » after pane- 

* 

Crating to low altitudes inside the cleft, could diffuse to adjacent field lines, 

the'.eby contributing to the polar rain. Figure Ic shows the case of a spectrum 
which consists of a power "lav and a Maxwellian population, but no Gaussian con- 

tribucion. 

The expression for the error, given in Appendix I, gives an upper bound of 
the mean square error of the logarithmic flxxx. In the vast majority of cases, 
this is of Che order of lO”*' or even less. The highest values encountered by 
us so far in a few examples are of the order of 3x10*^ . 


Considering this worst 
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CM*, ihft upper bound of the root seen square error is then 5.5 x 10*^ > This 
quantity would have to be compered with the sosllest values of the logsrithaie 
fluxes which are of the order of 4. Hence, we conclude that our snelytical fit 
is remarkably good. 


8 


XXX. implications of the Characteristic Paraaatars 

The axaaples report ad in this section are based on data fros the inverted«>V 
events observed by the Lov Energy Electron Analyzer (LEE) onboard the AE-D 
satellite. Ve have included 10s of data on either side of each invert ed-V event. 
The list of inverted-V events was conplled by C. S. Lin (private communication) 
based on .in analysis of the electron spectrograas. 

There vere tvo 1£E detectors on the AZ-S satellite, one at an angle of 7* 
vith respect to the satellite axis and the other at an angle of 60*. The angles 
of the tvo detectors vith respect to the geosagaetlc field varied cf course over 
an orbit. The flux component parallel to the geomagnetic field, however, was 
always the dominant one of the 7* detector, white the perpendicular component 
remained dominant for the flux measured by the 60* detector. Therefore, for 
purposes of identification, quantities relating to the 7* detector will be 
designated by the subscript|| and those relating to the 6C* detector by the 

a) Some Properties of the Power-Lav Populitioc 

One characteristic feature of the polar cap spectra is that in general 
they do not contain the power- law population at energies measured by the LET 
instrument, i.e., above 200eV Ic.f., Foster and Burrows, 1977], The reason 
is evidently the softness of the precipitating primary flux. Based on a simple 
model of a constant field-aligned upward directed electric field above the 
satellite Evans 11974] has shown that the low energy down-streaming power-lav 
population represents degraded backscatterec primary electrons and secondary 
electrons which have been reflected downward by this electric field. This 
interpretation has been confirmed by the more rigorous treatment of Staones 
[1978, 1981]. The latter has shown that for a precipitating primary beam with 
Gaussian energy spectrum peaked at 500 eV the pewer-lav spectrum lies below 
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•bout 120 eV which would nak« it invisible to the LEE instrunent. Bence, the 
absence of the |>owerolsw populstion from the LEE spectrus csn be censtdered s 
signature of soft prinsry precipitation characteristic of the polar rein. This 
is conflmed by Figures 2 and 3 in which the average peak energy of the Gaussian 
population and the average temperature of the t^xwellian with and without a 
power-law component are compared as functions of latitude at several local time 
sectors. The figures clearly show two interesting properties. The average peak 
energies and the average t«&peretures of the fluxes with power-law components 
are consistently higher than those without power-lav components. Secondly, both 
the peak energy and the temperature have a minimum In the latitude range between 
80* and 85* with the exception cf the temperatures in the noon (10.5 - 13.5 hrs) 
and afternoon (13.5 - 16.5 hrs) sectors. The slight increase in the mean value 
of T|| toward the top latitude range in those two time sectors is as yet unexplained. 

Table 1 lists the average and aaximun values of the Gaussian peak energies 
and of the temperatures separately for the fluxes with and without power-lav 
populations. Both the average and meximus values of the peak energies and of the 
temperatures of the 'fluxes without power-law populations are significantly lover 
than those of the fluxes with power-law populations. This is true of both the 
para21.el and perpendicular components. 

Figure 4 shows the average value of the power-law exponent c for eight 
local time sectors, for both the parallel and perpendicular fluxes. In most 
cases the magnitude of o has a caximua in the 80-85* latitude interval. 

Indicating that the rate at which the power-lav spectrum falls off generally 
increases toward higher latitudes. If the power-law population is sade up of 
secondary and backscattered electrons, this behavior of c is consistent with the 
increasing softness of the primary flux pclevard of the auroral zone. 
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b) Lost Cone Pictribtttlon of the Gausoion Beau 

Atnoaphara Explorer aatallitca have the capability either to apln about 

an axle which la aoatly perpendicular to the geonagnetlc field or to fly deapua 

aa detemlned by ground coBsaand [Spencer et al., 1973]. The eplanlng orblta 

provide an opportunity for an inveatigation of' the pitch angle variation of the 

varioua paranetera. Figurea 5 and 6 ahov a fev exas^les of the pitch angle 

varlatlona of the energy and the temperature T. These figures show 

aone Interesting features. Because of the invariance of the magnetic moment all 

Che dowBward hemispheres (0* ^ i ^ 90*) are filled with Gaussien fluxes at those 

relatively low altitudes. The upverd streaming fluxes have Gaussian popialstions 

in the pitch angle range from 90® to about 110* (in seme cases even to larger 

values) and no Gaussian fluxes beyond, k'hich is typical of loss cene distributions. 

Most of the doiaiward streaming Gaussian fluxes are very nearly isotropic. 

Because of the loss cone distribution of the upward streaming Gaussian beams they 

can be assumed to consist of mirrored particles. It would therefore be of special 

interest to know the apex angle of the less cone at the satellite altitude (and 

thus at all other altitudes). Unfortunately, this angle is difficiilt to determine 

from the data because the LEE spectrometer measures one spectrum per second, i.e.,^ 

# 

one spectrum only about every 25* of pitch angle. Bence, the data merely provide 
upper bounds of the apex angle (obtained from the largest pitch angle at which a 
Gaussian population has been observed) and lower bounds (obtained from the next 
larger pitch angle) . Table 2 lists the mirror altitudes corresponding to both of 
these angles from data of AE-D orbit 635.' The mirror points of most of the lower 
bounds lie below the surface of the earth and therefore provide no useful informa> 
tion. Ihe loT:€:^■t upper bound occurs at 127 km. Since the beam energy for this 
particular case is 2.35 keV, this altitude is close to where electrons of that 
energy are collisionally absorbed by the atmosphere [see e.g.. Banks et al., 1974]. 
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B«nee, th« ap«x anti* of tha loss cons of s 2.35 ksV bsaa is Indssd eloss to 
115* St ths sltltuds of 567 ^ (for an span angla of 114* tha nlrxoriag sltltuda 
is 162 te - too high for eollisionsl absorption of 2.35 kaV alaetrons). 

Figura 6 shows that tha taaparaturas ovar the downward hmisphsra (6 < 90*} 
and outside tha loss cone of the upward haslsp^ra (90 < 9 < 110 appr.) are 
slthar isotropic or vary only slightly with pitch angla. Inside tha loss cone 
the taoparsturas of the Maxwellian populations decrease (thus becosing highly 
anisotropic) with the minisnim occurring parallel to the field line. The ratio 
of T„ ^®***/T„ can be as large as 2 (and sonetimes even larger), 
c) Correlation Bettfeen the Gaussian Peak Energy and the leaperature T 

The energy at the ntaximua! of the Gaussian distribution appearing in equation 
(1) is given by 

E « 1/2 (E + ,'E 2 + )• 

p o o 

A correlation analysis of T and E has been perfonaend, and the results are 

P 

sunmariaed in Tables 3a and 3b. The correlation coefficients listed in Table 
3a have been calculated separatelj' for eech cell in magnetic local time and 
invariant latitude, while the coefficients listed in Table 3b have been cal> 
culated for each magnetic local time sector for all latitudes above 6C*. This 
breakdown is presented to study whether the correlation depends in some sig* 
nlf leant way on the magnetospheric region of origin of the fluxes. As the 
tables indicate, the correlation has a maximuc during daytime and falls off 
toward the night. The fluxes carried by the Gaussian populations behave in a 
very similar way, i.e., the average Gaussian fluxes (averaged over the lifetime 
of AE-S) have a maximue during daytime and decrease toward the night in all 
latitude intervals. Therefore, as the intensity of the Gaussian fluxes de- 
creases, their influence on the temperature declines compared with other mecha- 


nisms. 


I 
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The linear relation betveen T and. the peak energy. for all -.local tiaea 
and laticudof hao K«eo obtained a regression analysis vhicb yields 

• T|j • 0.604 + 0.400 E^ 

T, - 0.563 + 0.452 E^ . 

X 

The standard error of the regression coefficient is 0.0086 for the parallel 

coBponent and 0.0072 for the perpendicular one. 45% of the variance of T|| 

is explained by E • and 59% of the variance ez Ti is explained by E .* 

1^1 •*' 

Burch et el. [1976] have anticipated the correlation betveen and T based 

on an examination of a s&all sanple of spectra. Lin at; hoffcan [1979. 

perfomed regression analyses separately for each inverted-V event. They 

obtained a set of values for the slope centerec around a value of 0.2 td:ieh 

is half of our value. The reason for this discrepancy is believed to be due 

to the different definitions of T used by the two groups. . . 

Several years' ago Whalen and bicDiamid 119723 suggested that there nay be 

e distributed source of cold electrons throughout the acceleration region. 

The bean which has fallen through the entire potential drop would then act as 

a heat source of the lover energy population through collective effects [see 

for exanple, Shapiro, 19633, thus resulting in a correlation between the 

peak energy’ and the tet^erature of the Maxwellian population travelling parallel 

to the beam. Those Maxwellian particles travelling upward inside the loss cone, 

on the other hand, have probably been scattered up by the ambient electrons or 

ions and have thereby been cooled. Since the hiat transfer rate due to this 

sechamiss is proportional to the flux intensity of the beam, it is therefore 

also consistent with the reduced correlation between T and E^ observed during 

night tiae as discussed above. 
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d) Soae Properties of the One-Disensionel Velocity Dlstrlhutlon 

The shepe of the one-dlnenslonal electron velocity distribution is of 
special interest for plnsna physical reasons. In particular, if there la 

a huffip on Its tall, the plasma may be unstable. If the hump is being steadily 

% 

supplied with new particles, then the Instability can saturate at a non- 
equlllbrium quasi'eteady state of plasma turbulence which leads to anomalous 
transport effects. 

The existence of such effects has recently been deduced from stormtime 
data in the auroral none Ilontheia et al., 1976]. Another interesting consequence 
cf plasma turbulencct is the emission of a characteristic spectrum of electrostatic 
and electromagnetic waves. The one-dimensioztal distribution is defined as 







V ) . 
2 


where f(v , v , v ) is the three-dinensioael velocity distribution. For nanv 
X y 2 

applications the distribution of interest is the one-dimensional velocity distri- 
bution parallel to the geomagnetic field which is given by 

2t „ 

F(Vj, ) • / d?/ f(Vjj , Vj, , V) (2) 

o c 

where F is the azimuth in. the plane perpendicular to B^. The formal dependence 
of f(V||, ^ , F) on the azimuth 'i' is included for generality. In most real 
situations the distribution is expected to be cyllndrically symmetric about B^. 
Kaufmann and collaborators iKaufmann et al., 1976, 1978a, b; Kaufmasn, 1980] have 
.discussed the sienlflcance o.f the one-dimensional distribution fxxnction in great 
detail and hare also examined the compatibility of their observed distributions 
with various acceleration mechanisms. 
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Qafortunately, most InstruBcnts do not neasuro the three-diaenslonel velocity 
distribution. The usually observed quantity is the flux as a function of energy 
E in a solid angle elemext about a given direction, l.e., 

C>(E, 6, S') dE slnededV 

vhere 4»(E, 6, S') has the dimensions (length) (time) (energy) (solid angle)' 
and ^here 9 is the pitch angle. The desired three-dimensional velocity distri- 
bution f(v„, S') is related to the flux t(E, 6, S') by the expression 


f(v„. Vj^, « V dv„av e- - J it t 

^ vv./ + V * 


mv, 


«(E. ?, n 


v„2 + V 2 


dV|jd\^d? 
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vhere f(yjj > Vj_, S'l has the dimensions (leagth)”- G'elocityl"2 ^ jj, using this 
expression fcr f(vj| , S). in the calculation of the one-diaensional distribution 

F(V||) the flux C(E, 9, S) has to be understood in the sense that the parameters 
entering into 4'(E, 6, S'), as given by equation (1), must be considered to be 
functions of the pitch angle 6. 

The integration appearing in equation (2) is carried out numerically. Vhlle 
the lower limit of the integration over v^ can be taken as zero as required (since 
the integrand vanishes at » 0 according to ecrxation 3), the upper limit of 
the numerical integration is restricted to the largest v’aiue of for which a 
flux measurement exists for a fixed which is given by 

I’', I 
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where is the closest pitch angle to 90* ot which a spectrum was measured 

and is the maximum energy the instrument detects (20 keV in the ease of 

the LEE). 

If Cm/2) v,,^ < MAX (E , T. ) , then most of the contributions to the 
integral in equation (2), with the integrand given by (3) and CD» c^e from 
the energy range in the neighborhood of and T. Since the integrand 
decreases exponentially with E for E>MAX ^b^ * ottox introduced in 
the integral by leaving out the interval above is negligible provided 




If B V,. ^ /2 > MAX (E , T) , then the error is negligible provided 
I o 


B 


? I(v.)L + v,,2)>> |v,^ 


2 ‘'Xmax "M 


2 ’/i 


The numerical integration of equation (2) has to be carried out for each 

value of '’jj* The largest value of V|| for which T’Cv^j ) has been computed is 

6 X 10® cm/s which corresponds to an electron energy of 10 keV. The range 

0 < V, < (v.) is divided into 100 equidistant intervals resulting in a set 
— Jjr JJihBx ^ 

of V, . . For any fixed value of v.y the pitch angles corresponding to the set 

- Xk 

of V . are given by 6. « tan ^ , and the values of the seven parameters 

jjt ' k Vjj 

entering into the flux expression (1) are obtained by linear interpolation 
between their fitted values at the pitch angles at which the flux has been mea- 
sured. The value of the integrand is then calculated at each v^ ^ with the 
energy given by (m/2)(V||^ + integration appearing 

in equation (2) can be carried out with the integrand given by equation (3). 
This integration must of course be repeated for a series of values of to 
obtain r(vjj,) in the desired range of its arg\=ent. 
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lo order to get good coverege from e range of values of for the 
integrel In equation (2) the one-dimensional distribution function F(V|| ) 
must be calculated frrai spinning orbits by using data from both LEE 
instruments. Since there were only relatively few spinning orbits cf 
it was only possible to obtain a limited number of one-dimensional distri- 
bution functions. Some examples of different cases are presented in Figures 
7 through 10. Figure 7 shows a high velocity tail of a completely stable 
distribution. Figures 8 and 9 are two examples of plateau formation indicating 
the existence of an instability at a prior time, and Figure 10 shows a bump. 

An important special case e>:ists when the three-dimensional velocity distri- 
bution f(vj, j is spherically symmetric, i.e., if it depends on the velocity 

components vj^ and v^^ only through * Vj, ^ + Vj}. It is well known that' if a 
three-dimensional velocity distribution is spherically symmetric (or isotropic) , 
then the one-dimensional distribution as defined above can have at most one 
maximum and this is located at v « 0, irrespective of the shape of the three- 
dimensional distribution. Thus a hiamp or. the tail of the one-dimensional 
velocity distribution is automatically excluded if the three-dimensional distri- 
bution is spherically symmetric, even if the three-dimensional distribution has 
a hump as a function of energy. Although this is a well known property of 
distribution functions, a brief proof of this theorem will be presented in 
Appendix II because of the importance of this property for electron precipitation 
fluxes . 

As is indicated in Figure 5, the htimp in the energy distribution may be 
isotropic over the downward hemisphere. The above mentioned property of 
spherically sytssaetric three-dimensional distributions' applies separately to 
Vji^ 0 and ^ 0, and therefore the one-dimensional velocity distribution has 
no peak on the downw^ard tail. Since all spectra examined by us have picch angle 
distributions similar to that shown in Figure 5, it mey be generally true that 
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at low alcitudos the one-diaenslenal velocity dlccribucions of precipitating 
electrons generally do not have a huop on the high velocity tall. This le 
expected to be the case at altitudes below some critical altitude at which 
the beam particles «7lth the largest pitch angles first start to mirror. 
Bemlspherical Isotropy Is also the reason why Raufmann et al. (1978a;b] did 
not detect any secondary peaks » aside froo a few exceptions. 

As has been mentioned above, the existence of a hump on the tell of the 
one^dimenslonal distribution function may give rise to a plasma instability which, 
according to quasilinear plasma theory, has the effect of flattening the peak • 
into a plateau on the tail of the distribution. Papadopoulos and Coffey [1974a], 
however, have shown that under certain conditions a parametric instability 
is excited, the so-called oscillating two-stream instability. This instability 
has the effect of transferring wave energy out of the wave regime of phase 
velocities near the beam velocity (corresponding to frequencies near the plasma 
frequency) into the regime of much lower frequency ion fluctuations. This limits 
the growth of wave amplitudes of those waves with phase velocities near the beam 
velocity which strongly interact with the beam. As a result the beam is stabilized 
against quasilinear diffusion. According to Papadopoulos and Coffey J19J4aJ the 

condition for stability against quasi-llnear diffusion of the beam is 

2/3 7/3 2/3 1/3 . 1/3 

(n^/n^) (v^/Av^) ^ 


where m. and v. are the beam density and velocity respectively, n is the ambient 
electron density, iv.^ the velocity spread of the beam, the thermal velocity of 
the ambient electrons, M and m the ion and electron masses respectively, k the 

fit 

wavenumber with the ma>:inum growth rate for the parametric instability, and 

the Debye length. In none of the cases examined by us was this condition satisfied. 

As a result, the humped distribution function shown in Figure 10 is not expected 
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to stabilize into s qussi^stesdy ststo of plsssa turbulonc*. 

In those coses where the above condition is satisfied the growing ion 
density fluctuations give rise to an enhancement of the parallel resistivity 
of the plasma. Papadoupoulos and Coffey 11974b] have shown that in this ease 
the anomalous resistivity can be expressed in tetms of an effeetlye 
collision frequency 


n 

a 



eff 


(5) 


where is the electron plasma frequency of the ambient plasma. The effective 
collision frequency is given by Papadopoulos and Coffey [lS74b] 


V ‘ iv 

%£f - 

e e b 


( 6 ) 


where n, is the density of the. beam, n the densitv of the ambient electrons, 

V. the beam velocity, v the thermal velocitv of the ambient electrons, 4v. the 

0 6^ D 

velocity spread of the beam, and o is a factor of order unity. 

Once the one-dimensional distribution r(V||) is knotN'n, one can determine 
v^iether it has a maximum (i.e., a hump) on the tail, and, if so, fit another 
function of V|| to this hump, designated by T^Cv,,), representing the one- 
dimensional distribution of near mono-energetic beam particles. The density of 
the beam electrons is then given by n^ ■ / | ) , the beam velocity v^ 

is the velocity value at the maximum, and the velocity spread is simply the 
standard deviation of the distribution ). The ambient density n^ and the 

ambient temperature T^ (which determines the thermal velocity v^) are being 
measured by other instruments. Eence, all paratc-ters entering in equation (6) 
can either be measured or calculated, and the anomalous resistivity can thus be 
obtained for every spectrum which has a one-clsensional hump and satisfies 
condition (4). 
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XV. SuBBary 

A aathod for the eonputerized paraact trizaelon cf electron energy spectra 
has been outlined together with a few exanplcs of studies based on these derived 
parameters. It was shown that the presence of the pover*lav population is linked 

to the energy of the precipitating primary bean. This is consistent with the 

* 

% 

view that the power-law population is conpoaed' of secondaries and backseattered 
primaries. A detailed correlation study between and ? shoved that these two 
parameters are positively correlated and that the correlation is significantly 
stronger during local daytime than at nighttime. This effect may be related to 
the fact that the Gaussian fluxes also are larger during daytime and fall off 
toward the night. The one-dimensional velocity distributions parallel to the 
geomagnetic field are decreasing with in the great majority of cases, i.e. , 
they represent stable configurations. In some cases, however, they show a 
plateau on the tail, indicating an instability at an earlier time. A few hunp- 
on-the-tail distributions were also detected. 
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TABU 1: COMPABISOK OF TEMPEitATUBES AMD PEAK EKEHCIES IK THE 

PRESEKa AMD ABSEMCE OF POt^R-LAW FLUXES. 


^^ O )-0 

4 0 

- 0 

4 0 


Temperature (keV) Peak energy (keV) 


Average 

Maximum 

Average 

Maximum 

0.571 

11.978 

1.035 

3.637 

2.224 

19.990 

2.645 

11.950 

0.581 

4.705 

0.996 ' 

3.653 


2.344 


19.782 


2.753 


14.353 
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TABLE 2 

Upper and Lower Bounds of >iirror Heigh,.s of 
Electrons on the Surface of the Loss Cone 


UT 

Pitch * 
_ _ Angle 

Spacecraft 
Altitude (to) 

Mirror ** *** 
Altitude (to) 

Bean 

Energy (keV) 

Upper (U) or 
Lower (L) Bound 

13126 

110* 

529 

249 

4.58 

U 

127 

134* 

530 

0 

0 

L 

131 

• 

00 

M 

H 

534 

0 

0 

L 

132 

104* 

536 

399 

4.98 

U 

lAO 

92* 

544 

541 

5.68 

V 

141 

116' 

545 

71 

0 

L 

146 

121* 

550 

0 

0 

L 

147 

97* 

552 

518 

4.82 

V 

155 

99* 

560 

503 

1.73 

U 

156 

123* 

561 

0 

0 

L 

16: 

% 

139* 

566 

0 

0 

L 

161 

115* 

567 

127 

2.35 

V 

170 

105* 

577 

419 

4.61 

V 

171 

130* 

578 

0 

0 

L 

175 

132* 

582 

0 

0 

« 

L 

176 

108* 

583 

355 

4.72 

V 

185 

112* 

593 

251 

4.77 

V 

186 

137* 

594 

0 

0 

L 


* The pitch angle is taken with respect to the cevnverd direction. 

** If the cirror point would have fallen below the surface of the earth, 
it was recorded as 0 ks altitude. 

*** If no beas feature was present in the spectrur., the energy was recorded 
as 0 keV. 







TAB1.F. 3. CORREUTIOM CORFFICIENTS OP T AND R 
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Figure 1 


Figure 2 


Figure 3. 


Figure 4. 


Figure 5. 


Figure 6. 


a) Exasple of en analytical fit of a spectrua contistiag 
of the sum of a power^lav, a Maxvellian, and a Gaussian 
distribution. 

b) Example of an analytical fit of a spectrum consisting 
of the sum of a Maxwellian and a Gaussian distribution. 

c} Example of an analytical fit of a spectrum consisting 
of the Stan of a power-lav and a Maxwellian distribution. 
Variation of Gaussian peek energy with invariant latitude 
in eight local time sectors. 

Variation of temperature with invariant latitude In eight 
local time sectors. 

Variation of power-law exponent with invariant latitude in 
eight local time sectors. 

Pitch angle variation of the Gaussian peak energy (~^) • If a 

spectrum does not contain a near-monoenergetic peak (l.e. , a 

Gaussian population), E is listed as zero. The left-hand 

o 

ordinate scales refer to the solid curves and the right-hand 
ones to the dashed curves. The data were taken from AZ-D orbit 
no. 633 in the altitude range from 556 km to 605 Iol, invariant 
latitude range, from 76* to 73* and at MLT 20 hrs. 

Pitch angle variation of the temperature CT) • The left- 
hand ordinate scales refer to the solid curves and the right-hand 
ones to the dahsed curves. The data were taken from AZ-D orbit 
no. 635 and the same altitude, latitude and as the data of 
Figure 2. 



IFitun 7. 


Figure 8. 


Figure 9. 


Figure 10. 


Fleeae phyeieelly steble one^dieenslenel velocity 41etrlbutiea 
perallel ^e the geoeagnetie field derived froo obeerved electron 
energy epectre of AE-F orbit 635. 

One-diaensionel velocity distribution perallel to the geonagnetie 
field derived froB observed electron energy spectra of AE»D orbit 
647. The distribution presents e good example of plateau formation. 
One<-dlaensional velocity distribution parallel to the geomagnetic 
field derived from observed electron energy spectra of orbit 
647. The distribution presents a good example of plateau formation. 
One-^dimensional velocity distribution parallel to the geomagnetic 
fiei d derived from observed electros energy spectre of AI*C orbit 
635. The distribution shovs a clear example of a hump»cs-the-tail. 
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Appendix I. Curve Fitting Froeedure 

The curve fitting procedure involves c linear least squares fitting to 
the electron flux speetrua using a superposition of three functions. 




Where F is the electron flux and £ the energy. The fomts of the functions 
c 

used are: 

Fower law 


TciK) - ?! 3^^ * 


Maxwellian 


T AZ) - Cti. ^37? exp(- 

D 


where C is a constant, n is the density and is the temperature. 
Gaussian 2 

E - F, 

F^3(E) - I expI- ( “) + Fj] 


The actual fitting is performed in Icg^ space, using the functions in the 
forms: 

Fower lav 

- El + y-j log.Z 

where F^ ■ exp(Xj) 

Pj . Ej 

Ma>yellian 






The fitting is perforned in such e vay z’mt it is possible to fit a 
particular electron flux spectrus vith any one of the three given functions 
separately or vith any coobination of the three. The basic procedure is to 
co&pute the least squares fit, compute the error between the observed spectrum 
and computed spectrum and select the least squares fit with the smallest error. 
The error is given- by the formula, 

1 ^ 

ERROR • : {log,«lF^(Ej3 - log,.Ir^(E )])2 

K - h - 1 i-1 0 .L lu c X 

P 

where, K is the number of data points used in the fitting 

h is the number of non-zero parameters in the fit 
P 

is the computed fltmt 

F is the obseri'ed flux 
o 

E^ is the energy corresponding to the i^^ data point. 

The first step in the fitting routine is to fit the eitire spectrum to the 
power law function and the >!axwellian function separately and choose the b<Bst 
fit. Next the routine searches for the optimal fit for a superposition of these 
same two functions. This search is eccoaplished by fitting the data points 


i 


# 


eorrespooding te cncrglM through to the Mucuollian function, 

oubtrneting the results fron the observed flux values and fitting the rcaafnder 
to the power lav function. The two fits are then added, the error is conputed 
and eoiqi>ared to the previous least error and the best fit chosen. This procedure 
is repeated fer h ranging from k ■ 3 to k « 13. 

% 

% 

To check for the presence of a nonoenergetic peak (approximated by a Gaussian 
function) the spectrum representing the best fit from the above search is scanned 
in the energy range to *13 to locate three consecutive flux values which lie 
above the computed spectrum. If three or more such points are found, th^ are 
deleted from the observed spectrum and the procedure outlined above is repeated 
for this adjusted spectrum. The points lying above tre original best f: vhich 

were deleted are fit to a 'aussian function and added to the recraputed power lav 
plus liaxwellian fit. The error is then found and compared to the original best 
fit. The fit with the least error is then taken as the overall best fit to the 
data and the parameters are stored in the output file. 


Appendix II. One-Dincneionel Ruop of en Isotropic Thrce-DlBcnsionsl Distribution 
Let tbe one-dinensionsl distribution function F(v ) have s naxianm st v . 

Z ZM 

Let the three-diaensional distribution function be isotropic, f(v^), vhere 

+ P(v ) end f(v*) are related by the expression 

X y z z 


P(V^) ■ / 7 f(V*) 


(A.1) 


Theores: has slvays the value 0 irrespective of the shape of f(v^), and 

' 2X3 

this is the only value it can have. 

Proof ; 

A necessary condition for a naxinua of T(v) is 
^ z 


-r^ F(v ) 
ov* z' 
z 




cv 

V -CD z 

zm 


*v 


dv cv 
X y 


• 0 


zm 


2v ; / f <v2) ; 

“ -i i(v5) ^ 


dv dv ■ C 
X y 
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Let us ass\sBe that v ¥ 0, Then 

21S 


-« d(v^) 'V 


dv dv 
X y 


m 0 
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Introduce polar coordinates in the v^v^-plane 


V « V, cosS 
X * 


V “Vi smc 

y ^ 


v2 M y: 2 Ay 

zm 


4 


Then ve obtain (ainee 


■) 


d(v2) dtfj) 


• 2it . 

/ du / d6u fte* 4> V* ) ■ 0 

0 O a«i*) ^ “ 


/ 


d(^2) _d_ + vL) - 0 


d(^2) 


f(«) - f(v2 ) - 0 
Za 


Since f(«) ■ 0, this laplies 


f(v^ ) • 0. 

zm 


(A. 2) 


Because of the isotropy bf f it follows that Since 


zs 


the function f (v^) is by definition non -negative, it can vanish at the 


point v^ « only if 


f(v2) 

d(v^) ‘v^ 

v2 


(A.3) 


zat 


fience, we conclude that 


Z 


0 with V y 0 
zm 


Z2S 


only if 


f(v*)j • 0. 

d(v*) *v*« 

V* 

2XH 


In order to detenaine whether F(v ) actuelly has a maxisum at v we 

Z 7M 

have to look at the second derivative. 


d2 

dv 2 2 


• 2 — ^ / r — dv dv 

dv * i eV XV 

2 -* 2 


2» V ^ d(v.^) 

%d(V) 


»V> + V 


z 


•2t I£(v 2) + f(v 2)3 

z z cv z 

z 


• 2 t If (v 2 ) + 2 v 2 — 1 _ f (v 2 )] 


d(v2) 


which vanishes at v2 » v2 according to equations (A. 2) and (A.3). Bence, 

ZXD 

there cannot be a na^iinun at v 9 ^ C. If f(v_, v , v^) is sphericallv 

ZQ X y z 

symaetric and if F(v ) has a naxinun, it can occur only at v • 0 . q.e.d. 

z z 


